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Abstract 

Using the renormalization group motivated smoothing technique we study 
the semiclassical structure of the pure Yang- Mills vacuum. We carefully check 
that identified clusters of topological charge behave like instantons around 
their centers. Looking at distance distributions of clusters we find instanton- 
antiinstanton pair formation in the deconfined phase. We use suitably nor- 
malized gauge invariant field strength correlators to obtain first quantitative 
information on the relative color orientation of instantons. We give further 
evidence for a modified instanton profile in the confined phase and present a 
simple alternative model for it. 
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The last decade has witnessed tremendous success in the practical application of the 
instanton model to describe hadronic physics [|I|-£§. This model is a radical simplification 
of QCD, taking just one specific feature of Yang-Mills theory into account, the existence 
of instantons Q. Within this framework it is very natural to understand chiral symmetry 
breaking (and its restoration at the hadronic phase transition) and to solve the U^(l) 
problem 0, but confinement has defied an explanation so far. The model has been formu- 
lated to represent a practical tool for a wide range of non-perturbative calculations at zero 
and non-zero temperature. Instanton models have been developed ranging from the random 
(uncorrelated) ideal instanton liquid model RILM ||, used mainly for zero temperature, to 
the interacting (correlated) instanton liquid model IILM ||, in particular for T 7^ 0, when 
the temperature dependence of the interaction becomes important ||10|| . Even in the RILM 



there is some trace of interaction in so far as without some repulsive interaction [IT] the 
existence of a distinguished instanton size could not be postulated. 

Presently, some cornerstones of the IILM become subject of interest for lattice investiga- 
tions and get immediate phenomenological importance. In deep inelastic scattering (DIS), 
attractive instanton-antiinstanton (valley) configurations within a narrow range of relative 
orientations, and the small size flank of the size distribution d(p) are substantial for the event 
rate estimated for instanton mediated multiparticle processes Jl2| . Recently there was some 



effort to provide parts of this information, defining instanton size and distance distributions 



on the lattice, an enterprise which is under current debate |T3|-[T5 . 

The basic idea of instanton models is to reduce the infinite number of integration vari- 
ables in the path integral, to an integral over a manageable set of relatively few collective 
coordinates. Those are attributed to instantons forming an interacting system that is de- 
scribed by the grand-canonical ensemble with the partition function 

1 N+ r N - r 

Z= E WJ^U [dVtn(pt)]]l [dnjn(pj)]eM-Sint). (1) 

Here f2- = (zf,pl,U[) are the collective coordinates of the N + instantons (e* = +) and AL 
antiinstantons (e* = — ); z\ are the center positions, p\ the sizes and the SU(N C ) matrices 



the color orientations of (anti-)instantons. n(p) is the semiclassical instanton distribution 
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The instanton interaction Si nt can partly be taken into account in a cut-off single-instanton 
density function. The form 
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results from a mean-field approximation jT/]] due to the instanton-instanton and instanton- 
antiinstanton interaction, being repulsive on average, b is the one-loop coefficient of the QCD 
/^-function (6 = lliV c /3 for pure Yang-Mills theory), and v = (6 — 4)/2. This results in a 
packing fraction y p 2 " 2 = where the instanton interaction is measured by 7 2 = |^?r 2 
and (3 = [3(p) « 15 the average action per instanton. The color factor 7 2 represents the 
averaging over (relative) gauge orientations. 

To compute mesonic correlation functions in this framework, Bethe-Salpeter equa- 
tions [1^] can be formulated in the so-called random phase approximation. For this purpose 



it is assumed that the instanton orientation is random and the distribution of instanton 
centers is uncorrelated. Baryonic correlation functions are difficult to obtain in this way. 
For crude estimates of the condensates, the knowledge of the instanton size distribution d(p) 
(or simply the average p) and the distribution of distances (or the average nearest neigh- 
bor distance) is sufficient. Much of the knowledge about hadronic correlation functions |J 
and about the spectral density |19| of the Dirac operator, however, has been obtained nu- 
merically, from Monte Carlo simulations of the instanton-antiinstanton liquid. Within this 
method of calculation it is possible to incorporate information on pair distributions in Eu- 
clidean and color orientation space which originates from specific Ansatze or is provided by 
lattice observations. 

Very little was known, so far, about the relative color orientation. Concerning the quark 
mediated interactions, the fermionic determinant which enters a non- quenched instanton 



Monte Carlo || can be reasonably parametrized in terms of the zero-mode overlap matrix 
element Tja oc i{u re \Az)p + p~ / {2p + p~ + (Az) 2 ) 2 in order to sample intanton orientations 
a la Metropolis []. The gluonic part of the interaction is even less known, i.e. strongly 
dependent on the Ansatz of II or IA superposition. Usually it is taken care of by a sharp p- 
distribution (^) and by a hard-core cut-off Jll|| . The most important and remarkably Ansatz 
independent source of orientational interaction is the dipole-dipole interaction which acts 
exclusively between instanton and anti-instanton and does not contribute to the coefficient 
7 because the mean-field approximation neglects color correlations. 

In this letter we give first quantitative results on the color correlation in instanton- 
instanton and instanton-antiinstanton pairs in pure Yang-Mills theory, which exists entirely 
due to gluonic interaction. Further we present the pair distribution between topological 
clusters as a function of distance in 4D space. Finally, we critically examine the assumed 
color coherence within a topological cluster. The progress in this paper towards a more 
detailed description of the instanton vacuum structure would not be possible without the 
method of constrained smoothing (CS) which makes use of the concept of perfect lattice 



action |F2U| . It allows to study semiclassical configurations which can be far from being exact 
solutions of the Euclidean field equations. Our Monte Carlo configurations are created 
with the same action using the Metropolis method. The CS procedure, in contrast to 
cooling or smearing, is essentially a two-lattice procedure, where lattice configurations are 
related to each other by blocking (fine-to-coarse) and inverse blocking |20| (coarse-to-fine). 
The smoothed configuration on the fine lattice, emerging from inverse blocking a blocked 
configuration, can be considered as semiclassical, which embraces deformations from classical 
solutions due to classical and quantum interactions. In the present work we want to get 
information about these interactions. In particular, nearby instanton-antiinstanton (AI) 
pairs will not be annihilated by CS. The distance and color correlations between these 



1 u Te \ is the real 4-vector which represents the SU(2) matrix f7 re i and Az = z 
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objects (as well as between instanton-instanton pairs) is not known from rigorous analytical 
calculations. Another possible deviation that we expect is a non-classical behavior of the 
field strength (lack of coherence) within a topological cluster which could give rise to a 
modified profile of an instanton. In contrast to CS, cooling does not converge before locally 
classical fields have been obtained. Therefore, a critical assessment of instantons in the 
quantum ensemble seems to be impossible. After some number of iterations, cooling or 
smearing destroy instanton-antiinstanton pairs. Therefore, based on these algorithms, it has 
been attempted to obtain detailed information on pair distance and p-distributions from a 
backward extrapolation of the cooling/smoothing history [|Hjj . 

The second ingredient, which enables us to investigate the color correlation between 
clusters and to examine the coherence inside clusters, is borrowed from the stochastic vacuum 
approach pl| . In this approach, a gauge invariant field strength correlator (the non-local 
gluon condensate) has been proposed as a model independent characterization of the gauge 



theory vacuum. This correlation function has been studied on the lattice |22] using the 
cooling method. Presently we are calculating P31 the two-point correlator 



D^upa = ( Tr ( G^{xi) S(x 1 ,x 2 ) G pa (x 2 ) S(x 2 ,x 1 ) ) ) (4) 



where a Schwinger line phase factor S(xi,X2) = Pexp i Jf* dx^A^x 



must be inserted. 



We use CS instead of cooling and the clover construction of the field strength components 



G a . In a semiclassical approximation |p4* |, the leading contribution linear in d(p) contains 
just all expressions as given by the classical instanton field, for instance - in the regular 
gauge - G^ v replaced by 
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Here z and p have the meaning explained before, 77" is the 'tHooft symbol and R ab = 
1 r-jn, ^ T ajj T bjj+^ j g ^ e adjoin representation of the global color orientation U of the instan- 
ton. To obtain the correlator (Q), the instanton center z and radius p have to be averaged 
over with the density @ after the phase factor S(xi, x 2 ) has been evaluated. 



In the present work, once the centers of topological charge clusters are identified, the 
distribution of relative color orientations can be obtained by measuring ([|) between the cen- 
ters for pairs of clusters. The classical coherence of a cluster can be examined by measuring 
(||) between the center and some other point at a distance up to a few times p. While in the 
model the Schwinger lines for these cases are radial with respect to the instanton centers 
(and therefore S{x\,X2) = 1), the non-Abelian Schwinger phase cannot be neglected in the 
lattice configurations. 

In all computations presented here we have used a simplified SU (2) fixed point action on 



a 12 x 4 lattice |[25|| , for which CS as described in p6| can be carried out in a theoretically 



consistent manner. Such actions provide scale invariant ideal instantons 0] and suppress 



dislocations \2%. Instantons with p < 1.75 a (somewhat smaller than the blocked level lattice 
spacing) are instable and cannot be made visible by CS. Our simulations were carried out 
in the confinement phase {(3 = 1.4 corresponding to a temperature 0.71 T c ) and slightly 
above the critical temperature T c (/3 = 1.6, 1.15 T c ). To identify instantons in smoothed 
but true quantum configurations we first looked for clusters of topological charge formed by 
lattice points with a charge density of equal sign above a certain threshold g t hr = 0.015/a 4 at 
(3 = 1.4 and 0.005/a 4 at (3 = 1.6. Next, connected clusters are defined by a site-percolation 
routine. The search for local maxima g max = m&x x (q(x)) inside the clusters identifies the 
instanton centers 2j. The instanton size can be obtained from g max 



P* = ( 2 \\ ) ' ( 6 ) 

as suggested by the classical profile. For all clusters we have measured a cluster charge 
Qciustcr = J2xeciust C T and a cluster volume V cluster = Eseduster 1- Fi S- 1 ( to P) shows 

the histogram of instanton radii according to ([]). The left flank represents the probability 
to find high values of qWx- The cut-off at large p results from the chosen threshold value 
Qthr applied to the cluster search. The observed mean value in Fermi of pp=\A = 0.51 fm is 
in good accordance with [|I^] where a similar p definitions was applied. It is larger than p 



reported in ||14|| . In deconfinement we find pp=\.§ = 0.40 fm. Note that in the deconfinement 



an 0(4) asymmetric profile becomes important, which might explain the unexpected shape 
of the p distribution Q. Fig. f (bottom) shows how the observed clusters are clustering in the 
Vdustcr-Qciuster plane. The solid line represents the relation between the two cluster properties 
predicted according to the classical instanton profile (]5|). The chosen threshold topological 
density has defined the horizontal scale, the total (unit) charge per cluster determines the 
vertical one. The dashed line is the result of a (somewhat deliberately chosen) Gaussian 
profile of a unit topological charge object. This curve makes clear why no cluster charge 
Qciustcr > 0.4 could be localized inside an individual cluster. The different performance of 
the two profiles gives an indication that in the quantum vacuum the instanton shape might 



be modified similarly to the fremon model considered in Ref []17j . 

The coordinates of the cluster centers are used to measure histograms of relative distance 
for different pairings of cluster charges. The observed number of pairs per bin is compared 
to a random distance histogram. The random histogram was obtained by a simulation 
on 30000 configurations on the same lattice where the lattice coordinates of I (A) centers 
were chosen randomly all over the lattice with an appropriate density. The ratio is shown 
in Fig. 2 and represents the two-cluster distance distributions in both phases. (Here no 
distinction was made between time and spacelike separations). If there were no repulsion 
or attraction, C(d) would be equal to unity. C(d) > 1 means attractive interaction (pair 
formation), C(d) < 1 repulsion. The situation for confinement is shown in the upper part. 
For II (AA) pairs no deviation from a random distribution is observed for d > 3 a. For 
the IA pairs there is a somewhat stronger repulsion up to d = 3 a than for equal sign 
pairs. In the deconfinement (bottom) a peak is observed for the IA distance distribution, 
signaling the formation of IA pairs with a center-to-center distance of 2-3 lattice spacings 
(about 0.28 — 0.43 fm). For the same-sign case (II or AA) there is a slight repulsion at 



2 In present calculations we employ a characterization of cluster shape in the deconfinement that 
takes care of that. 
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short distances. This IA distance distribution might explain the drop of the topological 
susceptibility \ a ^ t ne deconfmement phase transition as resulting from a local screening of 
topological charge. In our case X/3=i.4 = 1-0 ■ 10~ 3 /a 4 dropped to X/3=i.6 = 5.8 • 10 _5 /a 4 over 
the transition. 

We have modified the gauge invariant field strength correlator to a normalized two- cluster 
overlap function for any pair of clusters 

_ (Tr (^(xi) S(xi,x 2 ) G^(x 2 ) S{x 2 ,xi))) p&thes 
(TriG^xt)*))-* (Tr (G TX (x 2 f))i 

were x\ and x 2 denote the positions of the cluster-maxima. The Schwinger lines were 
chosen randomly among the shortest pathes connecting Xi and x 2 and the average was 
done incoherently as indicated in (|7|). For a non-correlated ensemble of well-separated ideal 
instantons (II) or antiinstantons (AA) one would expect that O averages to zero. In fact, due 
to the normalization (|7|) takes values in the interval (— |, 1). This overlap is a measure for the 
normalized adjoint trace of the relative gauge orientation between the clusters, Tr ad j U re \ = 
(Tr U rc \) 2 — 1. Ojj/aa — ~ \ corresponding to Tr U Te i = 0, and O = 1 meaning U re \ = ±1 
(color parallelity) . Averaging with the normal Haar measure (no color correlation) would 
lead to (O) = 0. For an IA pair Oja is expected to be trivially zero, simply due to the 
opposite duality in rfrjff^ = 0, whatever the relative orientation is. In order to get access to 
the relative co lor orientation, in spite of this color blindness of (^), we perform a twist in the 
field strength of one of the clusters in the IA pair, changing E a — ► — E a before calculating 
the overlap (0). 

In order to define U re \, the relative orientation between the clusters, we introduce an 
SU(2) matrix U into (0) and try to find, for each pair of clusters, the maximum of the 
function Ojj : 

c ^ ( Tr ( G »A X l) U S (Zi,^) G I1V (X 2 ) S(X 2 ,X 1 ) t/ + )) p athes /o v 

t = maxC/rr = max r^- ? . 8 

u u (Tr(G pa ( Xl )*)) l 2 (Tr(G rA (x 2 ) 2 ))^ 

We find U as the rotational matrix needed to produce maximum overlap, i.e. to undo the 
relative color rotation U Te \. We found that O measured for a cluster pair can be parametrized 
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as O = | (4 cos 2 (a) — 1) in terms of cos(a) = | Tr U as found by maximization with very 
good accuracy, for instance ( O — | (4 cos 2 (a) — 1) ) = —0.02 ±0.14 for ft = 1.4. In about 
95% of all pairs of clusters, a rotation could be found turning F almost to unity. Only 
5.4% of all clusters in the confinement phase were so incoherent that they reached values 
of F < 0.8. This is indirect support for the picture that the vast majority of the identified 
topological clusters contain all field strength components highly coherent in their center. 
We show in Fig. 3 the distributions in | Tr U re \ for all pairs of clusters with a distance 
d < 4 a, found in the Monte Carlo sample in the confinement and in the deconfinement 
phase, compared with pairs with d < 5 a. The dashed line shows the distribution according 
to the normal Haar measure (no color correlation) which would lead to vanishing average 
of O. The result is surprising in so far as in the confinement phase, opposite charge pairs 
are nearly uncorrelated while equal charge pairs correlated towards color (anti-)parallelity. 
The orientational alignment becomes stronger in the deconfined phase and is there roughly 
equal for both types of pairs. We are presently trying to identify the interaction parameters 
of S'int in ([!]) by instanton Monte Carlo simulations. 

Finally we will examine the question to what extent the clusters themselves can be 
considered as classical, coherent fields also outside of their centers. It is assumed in the 
instanton model that, throughout the entire instanton, the color orientation of all field 
strength components follows the same U in (|5]). Corroborating what we have discussed in 
relation to Fig. 1 (bottom), we find here that the semiclassical picture has to be revised 
since the two-point field strength correlator inside a cluster has a profile which differs from 
the classical one. In order to define a correlation length within a topological cluster, we have 
modified the field strength correlator to a normalized cluster profile 

v _ Tr(G^(a:i) S(x 1 ,x 2 ) Gjjjh) S(x 2 ,x 1 )) 

Tr^^) 2 ) U 

pinning x\ to the center of a cluster and considering it as a function of distance d = \x 2 — Xi\. 
For these calculations we restricted ourselves to on-axis displacements x 2 — X\ such that no 
path averaging was necessary. No twist was applied to the field strength at x 2 such that 

9 



fields of opposite duality could not contribute. The resulting profile function is shown in Fig. 
4 in the confinement (left) and deconfinement phase (right). Before final conclusions can 
be drawn the statistics of analyzed clusters has to be accumulated. The dashed line drawn 
in both figures represents a fit by a Gaussian profile function V(x) = exp(—x 2 /2B). The 
squares show an average over the field strength profile of individual clusters described by 
(fciuster^) = exp(—x 2 /B)/(7r 2 B 2 ) , averaged over the empirical B distribution obtained from 
<Zmax- The solid line is the best fit to the average using the classical profile. This comparison 
substantiates our tentative model introduced to explain the Qciuster^Kiuster correlation in the 
confinement phase. Let us define a r.m.s. size (for weighting with respect to the topological 
density) of the instanton p rms = y/2 B. From the B distribution we obtain the average 
B = 2.71 a 2 at (3 = 1.4. If the string tension is taken to be \fa = 440 MeV we find at 
(3 = 1.4 a lattice spacing a = 0.23 fm which corresponds to an p rms = 0.536 fm. At (3 = 1.6 
we find an average B = 2.87 a 2 . An extrapolation using the 2-loop formula gives a = 0.142 
fm at (3 = 1.6 which results in p rms = 0.340 fm. These numbers are very similar to the 
average radius defined via p according to ([]) and its distribution, although the modification 
of the profile at large distances is of principal importance for the confinement phase. We 
have to stress the possibility of asymmetry effects changing the instantons in the deconfined 
phase. This was not yet taken into account in this feasibility study. We are currently 
studying shape and correlations by treating time and space directions separately. 

In conclusion we emphasize that the use of constrained smoothing in combination with 
the measurement of gauge invariant two-point field strength correlators opens the possibility 
to obtain essential input information for semi-classically motivated models of QCD which 
was not available before. 

We thank D. Diakonov and M. Miiller-Preussker for initializing discussions and continued 
interest. 
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FIGURES 

FIG. 1. Top: Instanton size distribution according to @ in lattice spacings. The threshold 
values were rfthr = 0.015 in the confinement (left) and 0.005 in the deconfinement phase (right) such 
that instantons with radius p > 2.523 a or 3.321 a are not visible. Bottom: average and variance 
of a scatterplot of Ociuster vs. Truster- The agreement with the curve predicted by the classical 
profile (solid line) is lost for V^i uster > 8. A Gaussian instanton profile (dashed line) fits the data 
better in the confinement phase. Note that only 6.4% of the total number of clusters analyzed had 
a cluster volume Fluster > 8. 

FIG. 2. Distance distribution C(d) for equal charge (II and A A) and opposite charge (AI) 
pairs in the confinement (top) and deconfinement phase (bottom). The peak in the deconfinement 
phase for AI pairs indicates pair formation. Results obtained from 8930 clusters (1000 independent 
configurations) at /3 = 1.4 and 8862 clusters (2000 configurations) at j3 = 1.6. 

FIG. 3. Distribution in Tr U ie \ for II(AA) and AI pairs in the confinement (top) and deconfine- 
ment phase (bottom) averaged over all relative displacements d < 4 a (fat line) and d < 5 a (thin 
line). The random distribution (dashed line) is shown for comparison. 

FIG. 4. Cluster profile V (circles) obtained according to @ from 100 independent configu- 
rations containing 0(1000) clusters at (3 = 1.4 (left) and containing O(300) clusters at j3 = 1.6 
(right). Solid lines are fits to the classical profile with free parameter p, dashed curves correspond 
to the best fit to a Gaussian profile with free parameter B. Squares represent the averages over 
Gaussian profiles according to a B distribution obtained from g max for each individual cluster. 
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